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Stochastic excitation of non-radial modes 

I. High-angular-degree p modes 
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^Observatoire de Paris, LESIA, CNRS UMR 8109, 92190 Meudon, France 

ABSTRACT 

Context. Turbulent motions in stellar convection zones generate acoustic energy, part ol which is then supplied to normal 
modes of the star. Their amplitudes result from a balance between the efficiencies of excitation and damping processes 
in the convection zones. 

Aims. We develop a formalism that provides the excitation rates of non-radial global modes excited by turbulent 
convection. As a first application, we estimate the impact of non-radial effects on excitation rates and amplitudes of 
high-angular-degree modes which are observed on the Sun. 

Methods. A model of stochastic excitation by turbulent convection has been developed to compute the excitation rates, 
and it has been successfully applied to solar radial modes (Samadi & Goupil 2001, Belkacem et al. 2006b). We generalize 
this approach to the case of non-radial global modes. This enables us to estimate the energy supplied to high-(^) acoustic 
modes. Qualitative arguments as well as numerical calculations are used to illustrate the results. 
Results. We find that non-radial effects for p modes are non-negligible; 

- for high-71 modes (i.e. typically n > 3) and for high values of £; the power supplied to the oscillations depends on the 
mode inertia. 

- for low-n modes, independent of the value of £, the excitation is dominated by the non-diagonal components of the 
Reynolds stress term. 

Conclusions. We carried out a numerical investigation of high-.^ p modes and we find that the validity of the present 
formalism is limited to ^ < 500 due to the spatial separation of scale assumption. Thus, a model for very high-£ p-mode 
excitation rates calls for further theoretical developments, however the formalism is valid for solar g modes, which will 
be investigated in a paper in preparation. 



such represents a mechanical source of excitation; the sec- 
ond is caused by the advection of turbulent fluctuations 
of entropy by turbulent motions (the so-called "entropy 
source ter m" ) and as it such represents a thermal source o f 
excitation (|Goldreich et al.llTOM : iStein fc NordlimdIIMl . 
Samadi & Goupil (2001, hereafter Paper I) proposed a gen- 
eralized formalism, taking the Reynolds and entropy fluc- 
tuation source terms into account. In this model, the source 
terms are vifritten as functions of the turbulent kinetic en- 
ergy spectrum and the temporal-correlation function. This 
allows us to investigate se veral possible models of turbu- 
lence (jSamadi et al.ll2003ai rbl). The results have been com- 
pared with GOLF data for radial modes, and the theo- 
retical values are found to be in g ood agreement with the 
observations l|Samadi et al.|[2003bl ). Part of the remaining 
discrepancies has been recently removed by taking into 
account the asymmetry introduced by turbulent plumes 
(Bclk accm et al, 2006a, b). 

In this paper we take an additional step in extending 
the ISamadi fc Goupil (j200lt ) formalism to the case of non- 
radial global modes. This will enable us to estimate the 
excitation rates for a wide variety of p and g modes ex- 
cited in different types of stars. The present model pro- 
vides the energy supplied to the modes by turbulence in 
inner, as well as outer, stellar convective regions, provided 
the turbulent model appropriate for the relevant region is 
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1. Introduction 

' Amplitudes of solar-like oscillations result from a balance 
. between stochastic excitation and damping in the outer- 
' most layers of the convection zone, which extends to near 
[ the surface of the star. Accurate measurements of the rate 
. at which acoustic energy is supplied to the solar p modes are 
' available from ground-based observations (GONG, BiSON) 
as well as from spacecraft (SOHO/GOLF and MDI). From 
those measurements and a comparison with theoretical 
models, it has been possible to demonstrate that excita- 
tion is due to eddy motions in the uppermost part of the 
convection zone and by advection of entropy fluctuations. 

Stochastic excitation of radial modes by turbulent 
convection has been investig ated by means of sev- 
eral semi-analytical approaches (iGoldreich fc KeelevI 1977t 



Goldreich et all 119941 : iBalmforthI Il992l : ISamadi & Gounill 



20011 ): they differ from each other in the nature of the as- 
sumed excitation sources, the assumed simplifications and 
approximations, and also by the way the turbulen t convec- 
tion is described (see reviews bv IStein et aDl2004l : iHoudekl 
120061) . Two major mechanisms have nevertheless been iden- 
tified to drive the resonant p modes of the stellar cavity: 
the first is related to the Reynolds stress tensor and as 
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used. Studies of stoch astic excitation of solar radial modes 
(jSamadi et al.|[2003ai rbl) have given us access mainly to the 
radial properties of turbulence. The present generalized for- 
malism enables us to take into account the horizontal prop- 
erties of turbulence (through the non-radial components of 
the Reynolds stress) in the outermost part of the convective 
zone. 



2. General formulation 

Following Paper I, we start from the perturbed momentum 
and continuity equation 



d{po + Pi)v 
dt 



V : [pavv) = 

-Vpi +pigo + Po9i +Pififi , (1) 



In the Sun, high-angular-degree v modes (as high as on e 
thousand) have been detected fe.g.. lKorzennik et al.ll2004) . 
From an observational point of view. lWoodard et al.l (|200l[ ) 
found that the energy supplied to the mode increases with 
but that above some high-^ value, w hich depend on the 
radial order n (see lWoodard et 311120011 Fig. 2), the energy 
decreases with increasing They mentioned the possibility 
of an unmodelled mechanism of damping. Hence one of the 
motivations of this work is to investigate such an issue. 
As a first step, we develop here a theoretical model of the 
stochastic excitation taking into account the ^-dependence 
of the source terms to seek a physical meaning for such a 
behaviour of the amplitudes. 



Modelling of the mechanisms responsible for the excita- 
tion of non-radial modes is not only useful for high-£ acous- 
tic modes but also for gravity modes, which arc intrinsically 
non-radial. As for p modes, g modes are stochastically ex- 
cited by turbulent convection; the main difference is that 
the dominant restoring force for g modes is buoyancy. We 
however stress that convective penetration i s another possi- 
ble e xcitation mechanism for g modes (e.g. iDintrans et al.l 
|2005|) . Such modes are trapped in the radiative interior of 
the Sun, so their detection promises to give a much better 
knowledge of the deep solar interior. However, such modes 
are evanescent in the convection zone; their amplitudes at 
the surface are very small and their detection remains con- 
troversial. A theoretical prediction of their amplitudes is 
thus an important issue. It requires an estimation of the ex- 
citation rates but also of the damping rates. Unlike p modes, 
the damping rates cannot be inferred from observations 
and this introduces considerable uncertainti es; e.g. theo- 
retical estimates o f the g-mode amplitudes ()Goughiri985t 
iKumar et al.lll996[ ) differ fr om each other by orders ofmag- 
nitudes, as pointed out bv IChristensen-DalsgaardI l|2002D . 
We thus stress that the present work focuses on the excita- 
tion rates - damping rates are not investigated. A specific 
study of gravity modes will be considered in a forthcoming 
paper. 



The paper is organized as follows: Sect. 2 introduces 
the general formalism, and a detailed derivation of the 
Reynolds and entropy source terms is pr ovided. In Sect. 3, 
we de monstrate that the formalism of ISamadi fc Goupill 
(|200lh is a special case and an asymptotic limit of the 
present model. In Sect. 4, we determine, using qualita- 
tive arguments, the different contributions to the excita- 
tion rates and identify the dominant terms involving the 
angular degree (£). Sect. 5 presents the numerical results: 
excitation rates are presented and discussed. Sect. 6 dis- 
cusses the limitations of the model and some conclusions 
are formulated in Sect. 7. 



^ + V.((po+Pi)^)=0 



(2) 



where p is the density, p is the pressure, and g is the grav- 
ity. The subscripts 1 and denotes Eulerian perturbations 
and equilibrium quantities respectively, except for velocity 
where the subscript 1 has been dropped for case of no- 
tation. In the following, the velocity field is split into two 
contributions, namely the oscillation velocity (vosc) and the 
turbulent velocity field (m), such that v = Vosc + u. For a 
given mode, the fluid displacement can be written as: 



SVa 



1 



(A(t)^(r)e-''"«*+c.c.) 



(3) 



where loq is the eigenfrequency, ^(t") is the displacement 
eigenfunction in the absence of turbulence, A{t) is the am- 
plitude due to the turbulent forcing, and c.c denotes the 
complex conjugate. 

The power ( P ) injected into the modes is related to the 
mean-squared amplitude ( < |Ap > ) by (see Paper I) 



P^ri< \A\^ > / w, 



' 



(4) 



where the operator <> denotes a statistical average per- 
formed on an infinite number of independent realizations, 
r] is the damping rate, and / is the mode inertia. 

We use the temporal WKB assumption, i.e. that A{t) 
is slowly varying with respect to the oscillation period, 
77 w d \nA{t)/dt <C uJo (see Paper I for details). 
Under this assumption, using Eq. ^ with Eq. ([T]) and 
Eq. ^ (see Paper I) yields: 



dAjt) 
dt 



(5) 



where d x is the volume clement and S — ~{ft + ^ht +gt) 
the excitation source terms. Temporal derivatives appear- 
ing in Eq. ([5]) are 

— the Reynolds stress contribution: 

^ = -|(V:(po««)), (6) 

where u is the turbulent component of the velocity field. 

— The entropy term: 



d_ 

dt 



as 



"df 



(7) 



where dst is the turbulent Lagrangian fluctuation of the 
entropy ( as — dpi /dst ) and pi denotes the Eulerian 
pressure fluctuations. 

This contribution represents the advection of entropy 
fluctuations by turbulent motion and as such is a 
thermal driving. Note that it has been shown by 
iBelkacem et al.l (j2006bl ) that this term is necessary to 
reproduce the maximum in the amplitude as a function 
of frequency in the case of solar radial p modes. 
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the fluctuating gravity term 

dgt dpigi 



dt 



dt 



(8) 



where gi is the fluctuation of gravity due to the tur- 
bulent field. This contribution can be shown to be neg- 
ligible and will not be considered in detail here for p 
modes. 

Several other excitation source terms appear in the right- 
hand side of Eq. ([T|). However, as shown in Paper I, their 
contributions arc negligible since they arc linear in terms 
of turbulent fluctuations^ 

From Eq. ([5]), one obtains the mean-squared amplitude 



<\A\'it)> 



d^rid^r2e 



dtidt2 



V{tl+t2)+ii^o{ti—t2) 



X ((r (r-i) • S{r,,h)) (C(r2) ■ S*{r,, h))) , (9) 

where the subscripts 1 and 2 denote two spatial and tem- 
poral locations. 

To proceed further, it is convenient to define the following 
coordinates: 



_ r2 + r-i 
2 

r2 - ri 



tl+t2 

2 



where Xq and to are the average space-time position and r 
and T are related to the local turbulence. 

In the following, Vq is the large-scale derivative asso- 
ciated with Xq, Vr is the small-scale one associated with 
r and the derivative operators Vi and V2 are associated 
with ri and r2 respectively. The mean-squared amplitude 
can be rewritten in terms of the new coordinates as 



(MP it)) = 



1 



X ( [C -Sixo - -,to 



(10) 



4(c^o/)2 

/•2(t-to) r 

X / dioe'''^*"-*) / dr / d^xo d^r e-''"°" 

h{to-t) J 

-]) {^■S*[^o + -^,to + - 

Subscripts 1 and 2 are the values taken at the spatial and 
temporal positions [xq — ^, —^] and [ajg-l-^, ^] respectively. 
In the excitation region, the eddy lifetime is much smaller 
than the oscillation lifetime 1/77) of p modes such that 
the integration over t can be extended to infinity. Hence 
all time integrations over r are understood to be performed 
over the range ] — 00, +oo[. 

We assume a stationary turbulence, therefore the source 
term (S) in Eq. (|10p is invariant to translation in tg. 
Integration over to in Eq. (jlOp and using the definition of S 
(Eq. (HD, Eq. and Eq. dHD) yields: 



< \A\^ >-- 



1 



{Ci + Cj + CRs) , (11) 



where C|, and C| are the turbulent Reynolds stress and 
entropy fluctuation contributions whose expressions are re- 
spectively: 



^ Linear terms are defined as the product of an equihbrium 
quantity and a fluctuating one. 



— the Reynolds source term: 



^2 



d'^.To 



dr e 



d> 



X (po u,u. Vie)' {Po uiu^ V[)r")") (12) 

where a separation of scales is assumed, i.e. that the 
spatial variation of the cigenfunctions is large compared 
to the typical length scale of turbulence (see Sect. [6] for 
a detailed discussion), 
the entropy contribution 



+ OC 

S - I d^xo / dre-'"^"^ I d^r 

— 00 



(13) 



Crs is the cross-source term, it represents interference be- 
tween source terms. For p modes, Crs turn out to be neg- 
ligible because it involves third-order correlation products 
whic h are small and strictly vanish under the QNA assump- 
tion (jBelkacem et al.ll2006bD . 



2.1. Turbulent Reynolds stress contribution 
Eq. (|12p is flrst rewritten as: 

f 4-00 



c 



R 



dre-"^"'' I d^rpl x 



(14) 



The fourth-order moment is then approxi mated assum- 
ing the quasi-normal approximation (QNA, (|Lesieuilll997L 
Chap. VII-2)) as in Paper I. The QNA is a convenient 
means to decompose the fourth-order velocity correlations 
in terms of a product of second-order velocity correlations; 
that is, one uses 

((WiWj)(l)(M(U,„)(2)) = {{u,Uj)(i)) {{uiU,n)(2)) 

+ ((Wt)(l)(Mi)(2)) ((Wi)(l)(M,„)(2)) 

+ ((w^)(l)(Um)(2)) ((Wj)(l)(w/)(2)) (15) 

A better approx i mation is the closure model with plumes 
(jBelkacem et al.l l2006al lbl) which can be adapted to the 
present formalism in order to take into account the presence 
of up and downdrafts in the solar convection zone. 

It then is possible to express the Fourier transform of 
the resulting second-order moments in term of the turbu- 
lent kinetic and entropy energy spectrum (see Paper I for 
details) 



4>ij = TF({u,Uj)) = . ■ \ Or 



47rfc2 



fc2 



(16) 



where E{k,uj) is the turbulent kinetic energy spectrum. 

The turbulent Reynolds term Eq. takes the follow- 
ing general expression under the assumption of isotropic 
turbulence: 



Go = TT 



d^xo {plb*jbim) S'(^™(wo) 



(17) 



4 

where 

nijlra 
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2.2. Entropy fluctuations contribution 



The entropy source term is computed as for the Reynolds 
contribution in Sect. 12.1] Then Eq. (fTH)) becomes 



' — oo 

X -■ Xk{UJO + Xk[UJ) 



fc4 

rpijhn / ^il _ 



k^k" 



Si ■ (Vo : ^) ■ Bj 



(18) 

(19) 
(20) 



(28) 



where 



B k^ 



with 



du)Xk{^^a +cj)xfe(w) 
(29) 



where {e,} are the spherical coordinate unit vectors, {k,iu) 
are the wavenumber and frequency of the turbulent eddies 
and turbulent kinetic energy spectrum E{k,uj), which is 
expres sed as the pr oduct E{k)xk{^) for isotropic turbu- , ^ /j \ • i / t\ i 

lence (lLesieudll997D . The kinetic energy spectrum E{k) is ^^crc Es{k) is the entropy spectrum (see Paper I), and 

normalized as 



dkE{k) = -<^uP' 



(21) 



h'^ = |Cr Vl(ln|a, |)V^2(ln|a3 |) 
-C*Vl(ln|a, \)Vi{C) 

-CVl(ln I a, \)V{{C*) + V\{C*)V{{C) , (30) 



where w is an estimate for the ver tical con vective velocity ^^'^^t ^ = Y" ^ ^^"t mode compressibility, 
and $ is a factor introduced by (|T97l to take The final expression for the contn^t-" 

into account anisotropy effects. A detailed discussion of t^ations reduces to (see Appendix ESI ) 



the temporal correlation function [xk] is addressed in 
iSamadi et all (l2003br ). 



The final expression for the contribution of entropy fluc- 
Lppend 

d^xoa^ (A, + Bf) 5s(c^o) (31) 



Ci 



The contribution of the Reynolds stress can thus be ^here Ti is the anisotropy factor introduced in Paper I 

which, in the current assumption (isotropic turbulence), is 
equal to 4/3. In addition. 



written as (see Appendix I A. 1 
dm R{r) Sr{ujo) 



(22) 



with 

R[r) 



16 




2 

H 


44 




2 




15 


dr 




" 15 


r 







r dr 



+ c.c 



2 

5^ 

r 



11 



i> |2 22^*^/1 
15 15 



dc a 



^55(^0) 

where 



d (In I cts 
din?' 



dDf_ 



d In r 



(32) 
(33) 



dk 

— E{k)E,{k) I dtoxk{LOo+Lo)xkH (34) 



dr r 
16 



+ c.c 



r dr 



(35) 



15 5 
where we have defined 

Cr ^ ^ + i(e.-a) 

dr r 



f,|m| 



(23) 3. The radial case 



(24) 
(25) 



e,\r 



SRiuJo) 



|to|(2£+ 1) 



{£{£ + 1) ^ {m^ + 1)) 



We show in this section that we recover the results of 
Paper I providing that: 

— we restrict ourselves to the radial case by setting £ = 

(a = 0). 

— we assume a plane-parallel atmosphere. 
In the entropy source term ( C| ) the mode compressibility 



dfc 



E^{k) duj Xkii^ + 1^0) Xkiuj) 



(26) for a radial mode becomes: 

5p ^ 1 d{r^^r) 



(27) 



C 



Ye 



Note that in the present work, nonradial effects are 
taken into account only through Eq. (j23p . A more complete 
description would require including anisotropic turbulence 
effects in Eq. (fT8|) . this is however beyond the scope of the 
present paper. 



and from Eq. (jA.45p and Eq. (|A.46|1 . one then has 

2 



Ae=o = 



1 



r^ 




Dr 



d 



din; 



■ In I as 



dDr 



dllli 



(36) 

(37) 
(38) 
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Wc thus obtain (Eq. (|A.44p ) 



C.I 



1 



din I a, I 



(39) 



d In r d In r 
For the Reynolds stress contribution, Eq. ((22)l reduces 



to 



Cl 



Ait" 



{- 




2 


44 




il5 


Ar 




^15 


r 



^f§^ + --)'5HM (40) 



To proceed further, we use the plane-parallel approxima- 
tion. It is justified (for p modes) by the fact that excitation 
takes place in the uppermost part of the convection zone 
( r/R « 1 ). It is valid when the condition r kog^ ^ 1 is ful- 
filled in the excitation region (fcosc being the local wavenum- 
ber), i.e. where excitation is dominant. Consequently: 



d^r 



dr 



> 



(41) 



Validity of this inequality has been numerically verified and 
is discussed in Sect. H (Eq. 

With Eq. gj), Eq. dall), and Eq. (gO]) simplify as: 



Ss{ojq) (42) 
(43) 



^2 4^'^ / ,3 2 

Cs = —2- / d xo a. 



Ci 



64 
15 



d^r d In I 



dr d In r 







j d^xo pI 


dr 



de. 

dlnr \^ dr 

2 



These are the exp ressions obtained by Paper I and 
ISamadi et al.l (j2005[ ) for the radial modes in a plane- 
parallel geometry. 



4. Horizontal effects on the Reynolds and entropy 
source terms 

Wc derive asymptotic expressions for the excitation source 
terms (Eq. and Eq. (|A.44[) ) in order to identify the ma- 
jor nonradial contributors to the excitation rates in the so- 
lar case. 



4.1. t dependence of the eigenfunctions 

P modes are mainly excited in a thin surface layer which is 
the super-adiabatic zone at the top of the convection zone. 

Let us consider the equation of continuity and the trans- 
verse component of the equation of motion for the oscilla- 
tions. Let us neglect the Lagrangian pressure variation and 
Eulerian gravitational potential variation at r = i? (the 
surface). The ratio of the horizontal to the vertical dis- 
placement a^_th£_smfaceboundary is then approximately 
given by (jUnno et al.lll989l p. 105): 



6. 



(44) 



where a is the dimensionless frequency defined by: 

9 Jl, 9 



GM 



UJ 



(45) 



where uj is the angular frequency of the mode, R the star 
radius and M its mass. 

Frequencies of solar p modes then range between a ~ 10 
and cr ~ 50 (z/ £ [1,5] mHz). Hence, for the solar oscilla- 
tions, one always has: 



(46) 

However Eq. and Eq. involve coefficients depend- 
ing on the angular degree {£). We then also consider the 
ratio 



(47) 



Equation. ((T7|) is of order of unity for i ^ a. For example, 
for a typical frequency of 3 mHz, one can not neglect the 



horizontal effect 





in front of 




T 




r 



for values of £ equal 



or larger than 30. 

In what follows, we introduce the complex number /, 
which is the degree of non-adiabaticity, defined by the re- 
lation: 

Note that / = 1 for adiabatic oscillations. 

Let now compare the derivatives. Under the same as- 
sumptions as above, neglecting the term in {p/ p)d{5p/p)/dr 
in the radial component of the equation of motion (stan- 
dard mechanical boundary condition), one gets near the 
surface: 



dr \ r 



(/ri)-i[a^ + 2+(L7a^-2)(/ri-l)] 

'^V(/ri). 



(49) 



Hence, we always have \d£,r/dr\ ^ l^r/rl in the excitation 
region (except near a node of d^r/dr). Similarly to Eq. (|44|) . 
one can assume: 



dr dr 



(50) 



In fact, comparing Eq. (j50p with the numerically-computed 
eigenfunctions shows that it holds even better than Eq. 
in the excitation region. 

Finally, we can group the different terms of Eq. and 
Eq. (|A.44p in four sets: 



^1 = 



^2 



d^r 



dr 





2 






r 






r 



(51) 
(52) 



d^r 



dr 

2 



>S'4 



Idle, 

r 



d^h 



dr 



d^h 



dr 



(53) 
(54) 
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The terms in Si are always negligible compared to the 
others. At fixed frequency ( a ) we have thus: 

S-i > 5*3 > 54 > S2 for £ < cr (55) 

5"! > 53 > 54 « S2 for £ « CT (56) 

Si ^ S3 ^ S2 > 5*4 for £ « CT^ (57) 

£"2 > 5*3 > S"! » 5*4 for e > cr^ (58) 

In conclusion, the contribution of the horizontal displace- 
ment terms (6*2, 5*3) begins to dominate the excitation for 
e » a2. 

4.2. Source terms as functions of I 

Reynolds stress contribution; 

We start by isolating non-radial effects in the range 
i G [0;500]. Note that the limit i = 500 is justified in 
Sect. 16.11 by the limit of validity of the present formalism. 
We investigate two cases, £ <C cr^ and t k. respectively. 



The condition for which 



a is satisfied for around the 



/ mode for ^ > 50 and in the gap between the g\ and / 
mode, for t < 50. 

Using the set of inequalities Eq. ([55)1 to Eq. for a 

typical frequency of 3 mHz (i.e. a « 30), R{r) (Eq.| 
becomes for high-n modes {£ ^ cr^): 



dr 



(59) 



Hence, for high-rt acoustic modes one can use Eq. ([55]) in- 
stead of Eq. P5|) , and in terms of the excitation source term, 
the formalism reduces to the radial case for £ < 500 and 
high-n modes. 

For low-n modes (^ « cr^, i.e. for instance cr w 10) some 
additional dependency must be retained (see Eq. ([57]) ). One 
gets 



R{r) 



16 
15 



16 

15' 



5 



dr r 



+ c.c 



(60) 



The additional terms correspond to the non-diagonal 
contributions of the tensor V : ^ appearing in the Reynolds 
stress term because we are in the range £ w cr^ (see 
Eq. ([57])). The radial and transverse components of the di- 
vergence of the displacement nearly cancel so that Sp/p 
takes its minimum values. This is due to the fact that they 
are nearly divergence free (minimum of dp/p + V ■ ^ = 0), 
i.e. 



d^ 
dr 



L 



2^h 







(61) 



As the divergence of the mode corresponds to the diagonal 
part of the tensor V : ^, one can then expect that the exci- 
tation rate decreases (through the terms in d^* /dr x^^f^/r in 
Eq. ([50)1 ). However, such a decrease is compensated by the 
non-radial component of the tensor (C^ in Eq. (|60p ). Thus, 
for low-71 p modes there is a balance between the effect 
of incompressibility that tends to diminish the efficiency 
of the excitation and the non-diagonal components of the 
tensor V : ^ which tend to increase it. 



Entropy contribution: 

Numerical investigations shows that the non-radial com- 
ponent of the entropy source term does not significantly 
affect the excitation rates except for £ > 1000, which is 
out of the domain of validity of the present formalism 
(see Sect. 16. 1|) . The non-radial effects appear through the 
mode compressibility, |-Dc|^ (Eq. (jA.46p ). From Eq. ([ST]) 
one can show that non-radial contributions play a non- 
negligeable role for low-n modes. However, such low- 
frequency modes are not sufficiently localized in the supera- 
diabatic zone, where the entropy source term is maximum, 
to be efficiently excited. 

5. Numerical estimations for a solar model 

5.1. Computation of the theoretical excitation rates 

In the following, we compute the excitation rates of p modes 
for a solar model. The rate (P) at which energy is injected 
into a mode per unit time is calculated according to the 
set of Eqs. pT|) - (fT3)l . The calculation thus reqiures the 
knowledge of four different types of quantities: 

1) Quantities that are related to the oscillation modes: the 
cigenfunctions (^r,6i) associated eigenfrequencies 

2) Quantities that are related to the spatial and time- 
averaged properties of the medium: the density ( po ) , 
the vertical velocity ( w ) , the entropy ( s ) , and as = 
dPi/ds. 

3) Quantities that contain information about spatial and 
temporal correlations of the convective fluctuations: 
E{k), Esik), and Xfe(w). 

4) A quantity that takes anisotropy into account: <I> mea- 
sures the a nisotropy of th e turbulence and is defined 
according to[GougB ([1973) (see also Paper I for details) 
as: 



$ = 



< > 



(62) 



where = 



l\ and Uh is the horizontal velocity. 



To be consistent with the current assumption of 
isotropic turbulence we assume $ = 3. 

Eigenfrequencies and cigenfunctions (in point 1) above) 
were comput e d usi ng the adiabatic pulsation code OSC 
(jBourv et all 119751 ). The solar structure model used 
for these computations (quantities in point 2) w as ob- 
taine d using the stellar evolu tion cod e CES AM (jMorell 
I1997I ) for the interior, and a iKurucd (|l993f ) model for 
the atmosphere. The interior-atmosphere match point 
was chosen at logr = 0.1 (above the convective 
envelope). The pulsation computations used the full 
model (interior-l- atmosphere). In the interior niodeh 
we used the OPAL opacities (jiglesias fc Rogers! I1996D 
extended to low temperat ures with the opacities of 
'Alexander & Ferguson' ('1994'), the CEFF equation of state 
(Christensen-Dals gaard fc Daeppenl |1992| ). Convection is 
included ac cording to a Bohm- Vitense mixing-length for- 
malism (see lSamadi et al1l2006L for details), from which w 
is computed. The 4> value is set to 2 in the calculation. 
This is not completely consistent as we assume isotropic 
turbulence (i.e. $ = 3). This does not however affect the 
conclusions of the present paper as all results on nonradial 
excitation rates are normalized to the radial ones. Note also 



K. Belkacem et al.: Stochastic excitation of non-radial modes. I. 



7 




butions, one can rewrite Eq. ^ as 



(mHz) 



Fig. 1. top: The rate (P) at which energy is suppHed 
to each £, n mode for £ = 50, 100, 300 is divided by 
the excitation rate ( Pradiai ) obtained for the £ = 0, n 
mode. Computation of the theoretical excitation rates 
are performed as explained in Sect. 15.11 bottom: Ratio 
P I /{P I)radiai whcrc / is thc modc inertia. 



that thc equilibrium model docs not include turbulent pres- 
sure. These two limitations are of small importance here as 
our investigation in this first work on nonradial modes re- 
mains essentially qualitative. 

Finally, for the quantities in point 3, the total kinetic 
energy contained in the turbulent kinetic spectrum {E{k)) 
is obtained following ISamadi et all (|2006f ). 



5.2. Excitation rates 

The rate ( P ) at which energy is supplied to the modes is 
plotted in Fig. [1] normalized to the radial excitation rate 
( Prad ) ■ It is seen that the higher the £, the more energy 
is supplied to the mode. This is explained by additional 
contributions (compared to the radial case) due to mode 
inertia, the spherical symmetry (departure from the plane- 
parallel assumption), and the contribution of horizontal ex- 
citation. Note that, as discussed in Sect. [3] (see Eq. (|4ip). 
the departure from the plane-parallel approximation is neg- 
ligible for p modes. Then, to discuss the other two contri- 



P 



81 



C^2 



\UR)\ 



(63) 



where |Cr(-R)| is taken at the photosphere. Note that both 
terms of the product (Eq. (|63p ) are independent of the 
normalization of the eigenfunctions. Thus, as shown by 
Eq. (|63p . thc power supplied to the modes is composed of 
two contributions which both depend on t. 

The first is due to the mode inertia, which is defined as 



/ = 



/ dm 1^1^ = / 



16-1 



L' 16. 



'podr. (64) 



High-£ modes present a lower inertia despite the con- 
tribution in Eq. (|64p because they are confined high in the 
atmosphere where the density is lower than in deeper layers. 

The second term of the product Eq. depends on 
the non-radial effects through the excitation source terms 
(Eq. (|A.44p and Eq. (US))). To investigate this quantity inde- 
pendent of the mode mass (defined as // \^r{R)f'), wc plot 
the ratio PI / {PI)radiai in Fig. [T] Excitation by horizontal 
turbulence plays a role at all values of I and in the same 
manner. One can then discuss two types of modes, namely 
low-n (< 3)and high-n (> 3) modes (see Fig. [T|). 

— For high-n modes, non-radial effects play a minor role 
in the excitation source terms. The dominant effect (see 
Fig. [1]) is due to the mode inertia as discussed above. 

— For lower values of n, there is a contribution to the 
excitation rates due to the horizontal terms in Eq. 



Thus, contrary to high-n modes, the term in R[r) 
(Eq. (|23p ) is no longer dominant in front of the terms involv- 
ing for low-order modes. Turbulence then supplies more 
energy to the low-frequency modes due to horizontal contri- 
butions, which explains the higher excitation rates for low-n 
modes as seen in Fig[T] We stress that there is still turbulent 
energy that is supplied to the modes despite their nearly 
divergence-free nature. For such modes the non-diagonal 
part of the tensor V : ^, which is related to thc shear of 
thc mode, compensates and dominates the diagonal part, 
which is related to the mode compression. 



5.3. Surface velocities 

Another quantity of interest is the theoretical surface ve- 
locity, which can be compared to observational data. We 
compute the mean-squared surface velocity for each modc 
according to the relation (Baudin et al.. ,2005[ ): 



(65) 



where M = I /^"^{li) is thc modc mass, li the height above 
the photosphere where oscillations are measured. T,j = rj/ir 
the modc linewidth at half maximum (in Hz), and v1 the 
mean square of the mode surface velocity. Equation (|65p 
involves the damping rates {Tu) inferred from observational 
data in the solar case for \ow-t modes (see I Baudin et al.l 
I2005L for details). We then assume that the damping rates 
are roughly thc same as for the £ = modes. Such an 
as sumption is supporte d for low-^ modes « 50) as shown 
bv lBarban et all ((2004D . 
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0-0001 I \ ^ ^ ^ ^ \ ^ ^ ^ ^ \ ^ ^ ^ ^ ^ ^ ^ ^ ^ ^ ^ ^ ^ ^ 

0.5 1.0 1.5 2-0 2.5 3.0 

I- (mHz) 



Fig. 2. Surface rms velocities of £ = 0, 20, 50 modes calcu- 
lated using Eq. (|65[) and normalized to the maximum veloc- 
ity of the radial modes (see t ext). Note that the damping 
rates are taken from GOLF (jBaudin et al.l[2005h and are 
chosen to be the same for all angular degrees {£). Three 
(7 error bars derived from GOLF are plotted on the £ ^ Q 
curve. 



Fig. [2] displays the surface velocities for £ ~ 0,20, and 
50. Note that the surface velocities arc normalized to the 
maximum velocity of the £ = modes (Vq ~ 8.5 cm s~^ 
using MLT). This choice is motivated by the fact that the 
absolute values of velocities depend on the convective model 
that is used and is certainly imperfect. However, its influ- 
ence disappears when considering differential effects. As an 
indication, 3 cr error bars esti mated from GOLF for the 
£ = modes are plotted (see iBaudin et al.l l2005l for de- 
tails). The differences between the radial and non-radial 
computations are indeed larger than the £ = uncertainties 
for £ > 20. For a more significant comparison, error bars for 
non-radial modes should be used but they are difficult to 
determine with confidence (work in progress). For £ larger 
than 50, we do not give surface velocities, as derived those 
here depend on the assumption of approximately constant 
damping rate that is not confirmed for £ > 50. 

When available, observational data should allow us to 
investigate the two regimes that have been emphasized in 
Sect. 15.21 namely the high- and low-n modes. 

6. Discussion 

6.1. The separation of scales 

The main assumption in this general formalism appears in 
Eq. PT|) . where it has been assumed that the spatial varia- 
tion of the eigenfunctions is large compared to the typical 
length scale of turbulence, leading to what we call the .sepa- 
ration of scales. In order to test this assumption, one must 
compare the oscillation wavelength to the turbulent one or 
equivalently the wav enumbers. To this end, we use the dis- 
persion relation (see lUnno et al.|[T989t ) 

fc,^ = ^^(l-^)(l-^) and kl^^ (66) 

where N is the buoyancy frequency, Se the Lamb frequency 
and kr,kfi the radial and horizontal oscillation wavenum- 
bers, respectively, and = £{£ +1). 



For the turbulent wavcnumbcr, we choose to use, as a lower 
limit; the convective wavcnumbcr fcconn = ^tt/Lc where Lc 
is the typical convective length scale. Thus, the assumption 
of separation of scales is fulfilled provided 

In Figs. [3] are plotted the ratios kr/kconv and kh/kconv re- 
spectively. Those plots focus on the uppermost part of the 
solar convection zone where most of the excitation takes 
place. The assumption of separation of scale is valid for 
the horizontal component of the oscillation since one has 
kh/kconv *C 1 (for £ < 500) in the region where excitation 
is dominant. However, we must recall that our criterion is 
ba sed on the mi x ing len gth to compute kconv As shown 
bv ISamadi et all ^OQM) using 3D numerical simulations, 
the convective length scale (computed using the CESAM 
code, see Sect. 15.11) must be multiplied by a factor around 
five in order to reproduce the injection scale (Lc) in the 
superadiabatic layers. Hence, for a more conservative crite- 
rion, we must then multiply the ratio kh/kconv by a factor 
of five, this leads to a ratio near unity for £ « 500 (see 
Fig. Thus, for higher values of the angular degree, the 
scparation-of-scale hypothesis becomes doubtful. 

Concerning the radial component of the oscillation 
wavcnumbcr, the limiting value of £ seems to be the same 
(i.e. £ = 500). Thus, we conclude that for modes of an- 
gular degree lower than 500 one can use the separation of 
scales assumption. For £ > 500 the characteristic length of 
the mode becomes smaller than the characteristic length 
Lc of the energy bearing eddies. Those modes will then be 
excited by turbulent eddies with length-scale smaller than 
Lc, i.e. lying in the turbulent cascade. These eddies inject 
less energy into the mode than the energy bearing eddies 
do, since they have less kinetic energy. We can then expect 
that - at fixed frequency - they received less energy from 
the turbulent eddies than the low-degree modes. A theoret- 
ical development is currently underway to properly treat 
the case of very high £ modes. 

6.2. The closure model 

A second approximation in the present formalism is the 
use of a closure model. The uppermost part of the convec- 
tion zone is a turbulent convective system composed of two 
flows (upward and downward) and the probability distribu- 
tion function of the fluctuations of the vertic al velocity an d 
temperature does not obey a Gaussian law (jLesieu 3 [1993). 
Thus, the use of the quasi-normal approximation (QNA, 
iMillionshchikovl (j 19411 )) which is exact for a normal distri- 
bution, is no longer rigorously corr ect. A more real i stic clo- 
sure model has been developed in iBelkacem et al.l ()2006a|) 
and can be easily be adapted for high-.^ modes. This al- 
ternative approach takes into account the existence of two 
flows (the up- and downdrafts) within the convection zone. 
However, the QNA is nevertheless often used for the sake 
of simplicity as is the case here. Note that, when using the 
closure model with plumes, it is no longer consistent to as- 
sume that the third-o rder velocity moments s trictly vanish, 
however as shown bv IBelkacem et al] ()2006al lbl) . their con- 
tribution is negligible in the sense that their effect is smaller 
than the accuracy of the presently available observational 
data. 
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1=300 
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Fig. 3. Top: Ratio of the horizontal oscillation wavenum- 
ber to the convective wavenumber {kh/kconv), versus the 
normalized radius {r/R). kconv is computed using the mix- 
ing length theory such that kconv ~ 27r/Lc (Lc is the mixing 
length) and k^ is computed using the dispersion relation 
Eg. (|66p . Note that the ratio kh/kconv is computed for a 
frequency of around v = i mHz, depending of the angular 
degree [1). Bottom: the same as in the top but for the 
ratio kr j kconv • 



6.3. Mode inertia 

We have shown that the excitation rates for high-£ and 
n modes are sensitive to the variation of the mode iner- 
tia (/). / depends on the structure of the stellar model 
and the properties of the ei genfunctions in these external 
regions. ISamadi et al.l ()2006D have shown that different lo- 
cal formulations of convection can change the mode inertia 
by a small amount. This sensitivity then affects the com- 
puted excitation rates (P). However, the changes induced 
in P are found to be smaller than the accuracy to which 
the mode excit ation rates are derived from the current ob - 
servations (see iBaudin et al.|[2005t iBelkacem et al.ir2006bl ) . 
Furthermore, concerning the way the modes are obtained, 
we have computed non-adiabatic eigenfunctions using the 
time-dependent formah sm of Gabriel for convection (see 
iGrigahcene et alll2005l ). The mode inertia obtained with 
these non-adiabatic eigenfunctions exhibits a v dependency 
different from those obtained using adiabatic eigenfunctions 
(the approximation adopted in the present paper). On the 
other hand, the mode inertia using non-adiabatic eigen- 



fmictions fsee lHoudek et aLlll999l for details) obtained ac- 
cording to Cough's time-dependent formalism of convection 
(^Gough 1977) shows smaller differences with the adiabatic 
mode inertia. Accordingly, the way the interaction of os- 
cillation and time-dependent convection is modelled affects 
the eigenfunctions differently. As explained in Sect. 15.31 the 
formalism developed in this paper can be an efficient tool to 
derive constraints on the mode inertia to discriminate be- 
tween the different treatments of convection. Further work 
is thus needed on that issue. 

7. Conclusions 



We have extended the lSamadi fc Gout)ill ()200lD formalism 
in order to predict the amount of energy that is supplied to 
non-radial modes. In this paper, we focused on high-£ acous- 
tic modes with a particular emphasis on the solar case. The 
validity of the present formalism is limited to values of the 
angular degree lower than £ = 500, due to the separation of 
scale assumption that is discussed above, in Sect. 16.11 We 
have demonstrated that non-radial effects are due to two 
contributions, namely the effect of inertia that prevails for 
high-order modes (n > 3) and non-radial contributions in 
the Reynolds source term in (see Eq. (|22l) ') that domi- 
nate over the r adial one for low - order m odes (n < 3). 

Contrary to lBelkacem et akl ()2006b[ ) who used 3D simu- 
lations to build an equilibrium model, we restrict ourselves 
to the use of a simple classical ID MLT equilibrium model. 
Indeed, we were interested in this work in deriving quali- 
tative conclusions on nonradial contributions. Forthcoming 
quantitative studies will have to use more realistic equi- 
librium models, particularly for the convection description, 
such as models including turbulent pressure fe.g.^Bal mforthl 
|1992[ ) or patched models (e.g. [Rosenthal e t al. 1999). 

From a theoretical point of view, several improvements 
and extensions of the present formalism remain to be car- 
ried out. For instance, one must relax the assumption of the 
separation of scales if one wants to use very high-£ modes. 
Such an investigation (which is currently underway) should 
enable us to draw conclusions about the observational evi- 
dence that beyond some large value of £ the energy supplied 
to th e modes decreases with frequency (see IWoodard et al.l 
I2OOII . Fig. 2). Another hypothesis is the isotropic turbu- 
lence that has been assumed in the present work as a first 
approximation. Such an assumption is to be given up to 
get a better description of nonradial excitation of modes 
by turbulent convection, which requires further theoretical 
developments (work in progress). 

The present work focuses onp modes, but the formalism 
is valid for both p and g modes. We will address the analysis 
of gravity modes in a forthcoming paper. 
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Appendix A: Detailed expressions for source terms 

The eigenfunctions {$) are developed in spherical coor- 
dinates {er,eg,e^) and expanded in spherical harmonics. 
Hence the fluid displacement eigenfunction for a mode with 
given n, £, m is written as 



with 



i{r) = (ere,. + Vh) Yi,. 
d 



' H 



89 



1 d 



(A.l) 



(A.2) 



smf 



where the spherical harmonics {Yi^miO, 4>)) are normalized 
according to 

^y,,™r;„ = i (A.3) 

with ^l being the solid angle (dil = sin9d6d<j>) . 

The large-scale gradient Vq appearing in Eq. (fT3|) and p4|) 

for instance is given, in the local spherical coordinates, by 



d 1„ 

or r 



(A.4) 



A.l. Contribution of the turbulent Reynolds stress 



The Reynolds stress contribution can be written as (see 
Sect.O) 



C 



H - n I d^XQ [plb*jbim) I d'^k I duj 



E^k) 
~k^ 



k'k^ 



fc2 



where 

rpijlm _ ^il 

and 

bij = e, ■ (Vo : i) ■ 

where the double dot denotes the tensor product. 

We now consider the covariant (a*", a*, a"^) and the con- 
travariant (ar.,ae,a0) natural base coordinates where the 
eigenfunction can be expanded 



(A.5) 

(A.6) 
(A.7) 



(A., 



The natural and physical coordinates are related to each 
other by 

Cj = — ai (A. 9) 



where gij is th e metric tensor in spherical coordinates (see 
Table 6.5-1 in iKorn fc Kor nl lT96M ). i.e. 



1 



, = r'^ sin^ 9 , 



ffy = Ofori . (A. 10) 

Equation (|A.7p . with the help of Eq. (jA.Sp . can then be 
developed in covariant coordinates 

9^ _ „7 f 



da 



fdar\ 
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where T^^ is the ChristofFel three-index symbol of the sec- where 



ond kind l|Korn fc Kornlll968[) . According to Eq. (|A.lip and 
Eq. (|A.9|) . hj (Eq. (|A.7p ) can be written as: 



(A.12) 



To proceed, one has to express Eq. (|A.12p in terms of the 
physical coordinates (^'^). With the help of Eq. (jA.8j) and 
Eq. (|A.9[) . we relate the covariant coordinates qi to the phys- 
ical (contravariant) components 



(A.13) 



where the component ^'^ are derived from Eq. (jA.l[) 



dY, 



de 



1 ay. 



Hence Eq. (jA.12[) becomes 
/da\ dYi 

brd> ~ 



dr 



dr 



961 

1 dY,,^ 



her = 



Yo 



1 



r 

1 dY 



sinfl 
1 dYf 



r sm 6' 90 





1 


/92y,,„\ 


cos 6* 




r 


.sin2 6' 


902 j 


smc/ 





i?2 



i?, = 2 



df] /■ df7 



47r / 47r 



dr? /■ df7 



47r / 47r 



id 



kjki 
"fc2" 

^2 ^/ 



(A.18) 



We assume isotropic turbulence, hence the k compo- 
nents satisfy 



d^k = 



with Sij is the Kronecker symbol for i,j = r, 9, 0. Then we 
(A- 14) obtain 



i?2 



i?3 



^/l? (i: 



/-dr! fdl^fc /C2 ^ 2 



2(3 



(E|5.P + E(^s^.. + -) 



13 Ri + 2 (3 



where we have set 

= /"^ ^ R= M 

" ~ 7 4^ k'^ ' ^~ J 47r fc4 
To compute Ri , i?2 and i?3 , we write 
Brr = brr', Bee = bee', Baja, = 6aa 



(A.19) 



(A.20) 



Br 



1 ^ 9y,,,„ 

2 90 

2 '"'sin 6* 90 



The contribution of the Reynolds stress can thus be written j^^^ — — 
as: 



C| = 47r3y" dm J dk J duj R{r,k) 
E'^{k) 



with 



dr r 



(A.21) 



(A.22) 



fc2 



-Xfe(w + Wo)xA;(w) , 



Using the expression Eq. (jA.2ip for the quantities Bij , 
(A. 15) we obtain after some manipulation: 



where we have defined dm = Anr^podr and 



Ri 



Rir.k)^ I- 



dn f d^k 

4 "ij "im 



de. 


2 

+ 4 


k 


dr 




r 



2L\L^ - 1) 



Im 



rplj 



(A.16) 



+ \Crr~2i^+c.c) 



(A.23) 



Because T'J'™ = T-?™', it is easy to show that 

R{r,k) ^1^1 ^B:^B,,^ ^T^,lra^r^^,ml^ 

where B^j = (l/2)(6y -I- h^i). 

Using the expression Eq. (|A.6p for T*-''™, we write 



For i?3, some lengthy manipulation leads to: 

2 



d^r 


2 

+ 8 


k. 


dr 




r 



R{r, k) = R1-R2 + R3 



(A.17) 



1^2 -2(^^ + 2% + C.C) 
\ dr r r^ 

r dr 



(A.24) 
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where we have defined 
(2£+l) 
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A.2. Contribution of entropy fluctuations 



A 



ml + l)Ai 



1 



s We start from Eq. (|28p . and to proceed further in the deriva- 

l)-^£,|m|(p-25) tion of tire entropy fluctuation source term one has to com- 
pute 



Atn = ^W^+l)-m(m + l)) 



(A.26) 
(A.27) 



(A.38) 



To derive Ri, R2, and the following relations have been 
used 



Then, ^ and k are expanded in spherical coorindates 
{ar,ag,a^). We assume an isotropic turbulence, as a con- 
sequence the quantities kg, k^k^, kgk^ vanish when inte- 
grated over rife. One next obtains 



. COS0 dYe, 



96*2 sin 61 09 sin2 6' dcj)^ 

COS0, 



-L'-Yi^m (A.28) 



dfik h'^ = 27rH (Kr + hgg + h^^) 



(A.39) 



sinO 
dYtrn 



Ye^rn = ^£,mi^£,m+ie ^'^ + Be^rnYe,m-ie^'^ (A. 29) equal to 4/3. Assuming that as = cis{r) we have according 



where Ti. is the anisotropy factor introduced in Paper I 
which in the current assumption (isotropic turbulence) is 
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Ai,mYe,m+ie-"^ - i3f,™i^f,™-ie'^ (A.30) 



^ (Vffr,:„ • VHY,n) - 



in 
dfl 
An 



\sin9 d(j) 



(A.31) 
(A.32) 

(A.33) 



Combi ning Eq. (fOSl) . Eq. (fOO)) . Eq. (|X2l), and 
Eq. (|A.17p . with a = 1/3,/? =1/5 (isotropic turbulence, 
see Paper I for details), yields 
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For radial modes, this reduces to 
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(A.34) 



(A.35) 



The final expression for the Reynolds source term is 
then given by: 

Cl = 4^3 y a™ R{r) Sr^uj^) , (A.36) 
with 

/dk f 
E^{k) du xkito + uJo)xk{uj) (A.37) 

and R{r) by Eq. ((X34)) . 



to Eq. ([XI]) 
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(A.40) 



r^ sin^ ( 



8C 



To proceed, it is necessary to express the divergence of 
the eigenfunction 

C = Vo- i = DtY[^ (A.41) 
with 

D,{r, £) = Dr-—^h ; Dr = \^ {r^^r) , (A.42) 
r r'^ (jr 

where again = £{£ +1). 

We next integrate Eq. over dri/47r, the solid angle as- 
sociated with the eigenfunctions ^. One obtains with the 
help of Eq. (|A.4ip and according to Eq. (|A.40|) 

— / dfik n, = 



47r 
2nn 



\Df 
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din I Ug 
dlnr 



dDf 



dlui 



(A.43) 



The final expression for the contribution of entropy fluc- 
tuations reduces to: 

.2 _ 4^3 ^ [ ^3^^ ^2 ^ ^^^^^^ ^^^^^^ 



where 7i is the anisotropy factor introduced in Paper I that 
in the current assumption (isotropic turbulence) is equal to 
4/3. In addition, 

d(ln I P -"^ ^ 
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SsioJo) 
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(A.45) 
(A.46) 

(A.47) 



